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The spin 1/2 XXZ chain driven out of equilibrium by coupling with boundary reservoirs targeting 
perpendicular spin orientations in XY plane, is investigated. The existence of an anomaly in the 
nonequilibrium steady state (NESS) at the isotropic point A = 1 is demonstrated both in the weak 
coupling and weak driving limits. The nature of the anomaly is studied analytically by calculating 
exact NESS for small system sizes, and investigating steady currents. The spin current at the points 
A = ±1 has a singularity which leads to a current discontinuity when either driving or coupling 
vanish, and the current of energy develops a twin peak anomaly. The character of singularity is 
shown to depend qualitatively on whether the system size is even or odd. 

PACS numbers: 

I. INTRODUCTION 

The coupling of a quantum system with an environment or under a continuous measurement is often modelled 
via a quantum Master equation in the Lindblad form [l|, Q . Unlike a quantum system evolving coherently, the time 
evolution of which depends on the initial state, a quantum system with pumping tends to a nonequilibrium steady 
state, independent of the initial conditions. Recently there has been a large interest in investigating chains of quantum 
spins driven towards nonequilibrium steady states by applying gradients (pumping) at the edges. An important role 
in these studies is played by the one dimensional driven XXZ chain. The equilibrium XXZ model remains a subject 
of intensive investigations via various methods |^-|8|. In experiments, many transport characteristics in quasi one- 
dimensional spin chain materials can be measured Q. In the context of dissipative dynamics, recent numerical and 
analytic studies of a driven XXZ chain show a number of remarkable features like long range order and quantum 
phase transitions far from equilibrium [loj . negative differential conductivity of spin and energy 
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Recently it was shown that an isotropic point A = 1, a weakly driven XXZ spin chain has an anomalous spin 
transport [l6| . The current-carrying nonequilibrium steady state was obtained by bringing an XXZ spin chain in 
contact with reservoirs aligning the boundary spins along the Z-axis. The boundary gradient was small, which allows 
to think in terms of linear response of a system to the weak perturbation along the anisotropy axis. 



From a standard theory of linear response |17j one expects that the result of linear response will not depend on 
a form of boundary perturbation as long as it remains sufficiently weak. Respectively, an anomaly seen with one 
type of perturbation, should be visible with another type of perturbation as well. On the other hand, the linear 
response theory is valid for system in thermodynamic equilibrium at a certain temperature, while the driven chains 
with dissipative boundary driving are intrinsically nonequilibrium systems. At "zero gradient" the NESS describes 
a system without currents, in a maximally mixed quantum state, which corresponds to a Gibbs ensemble with an 
infinite temperature. To test further the linear response regime in a situation of a nonequilibrium, we investigate 
a weakly driven XXZ chain, where the weak boundary gradient at the boundaries is applied in the XY direction, 
transverse to the anisotropy axis m, 01 ■ Note that there are two complementary ways of realizing a weak drive in a 
Master equation formalism: (i) to set a weak boundary gradient, but finite or strong coupling to boundary reservoirs 
(ii) to set a weak coupling to reservoir, but finite or strong boundary gradient. 

Interestingly, in both contexts (cither weak gradient k <C 1 or weak coupling T <C 1) we see a distinct anomaly in 
the NESS. The anomaly becomes a singularity as T — > or k — > 0. of the magnetization current in the isotropic point 
A = 1. For weak coupling the anomaly sets in for N > 2 and for weak driving for N > 4. The presence of an anomaly 
is seen in practically all observables. In present communication we focus on singular behaviour of the steady currents 
of spin j and energy J E , and one-point correlations at the boundaries. The full analytic treatment is reported for 
spin current. In the limit k — > or T — > we observe a non-analyticity in the dependence J (A): the spin current is 
an analytic smooth function of A everywhere except at isotropic point A = 1 where its value is different. The nature 
of singularity, in addition, depends on the parity of system size, namely on whether the system size is an even or an 
odd number. 

Other observables e.g. current of ene rgy , density profiles etc. also have a singular behaviour at A = ±1. Thus, we 
can complement the result of Znidaric [ 1 61 ] with prediction of further anomalies at the isotropic point. In particular, 
we demonstrate that the energy current also develops an anomaly at the isotropic point, but with a different scaling 
and of a qualitatively different form. 
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The plan of the paper is the following. We introduce the model in the Sec{TlJ The cases of weak coupling and weak 
driving are treated separately in Sec lIIII and Sec lIVI In the conclusion we summarize our findings. 

II. THE MODEL 

We study an open chain of N quantum spins in contact with boundary reservoirs. The time evolution of the reduced 
density matrix p is described by a quantum Master equation in the Lindblad form Q (here and below we set 

h = lj 

^ = -i[H,p} + r L £ L [p} + r R £ R [ P }, (i) 

where H is the Hamiltonian of an open XX Z spin chain with an anisotropy A 

JV-l 

H = E (°t°t+i + °M+i + A^4+i) , (2) 
fc=i 

while C l and Cr are Lindblad dissipators favoring a relaxation of boundary spins k = 1 and k = N towards states 
described by density matrices acting on a single site p^ and pr, i.e. Cl\pl] = and Cr[pr] = 0. The Lindblad 
actions Cl and Cr are chosen in the form 

£r\p\ = "i YL=1 {P, VXV m } + ELl V m pVX, 

where W m and V m are polarization targeting Lindblad operators, which act on the first and on the last spin respectively, 

Wi = \/^F(of +iaf),W 2 = t/^*(o? -iaf),V 1 = J^(a v N + ia z N ),V 2 = - iafr). In absence of the 

unitary term in (pj the boundary spins relax with a characteristic time T^r^ 1 [23| to specific states described via 
the one-site density matrices pl and pr, satisfying Cl[pl] = and Cr\pr] = 0, where 

PL = \ (I - ko\) 
PR, = \ {I + k<t%) . 

From the definition of a one-site density matrix p on e-site = \ (J + XX '") '")' we see that the Lindblad superoperators 
Cl and Cr indeed try to impose a twisting angle of 7r/2 in XY plane between the first and the last spin j24| . The 
twisting gradient drives the system in a steady-state with currents. In the following we restrict to a symmetric choice 

r L = T R = r (3) 



for spin current while ^ Fr is chosen for studying the energy current 25]. The parameter k determines the 
amplitude of the gradient between the left and right boundary, and therefore plays a fundamentally different role 
from the coupling strength T. The limits k = 1 and k <C 1 will be referred to as the strong driving and weak driving 
case respectively. The two limits describe different physical situations as exemplified in other NESS studies. For 
symmetric choice ([3]) the steady state has a global symmetry [l!|,[22| 

p(N, A) = H x U rot Rp(N, A)RU+ t Z x , (4) 

where R(A <S> B <g> ... ® C) = (C ® .... ® B ® A)R is a left-right reflection, and the diagonal matrix U ro t = diag(l, 

is a rotation in XY plane, U rot (J%U^ ot = af l7 Urota^U^ = —cr^, and = (a x )® N . In addition, there are further 

transformations, see [l9| for details, relating the NESS for positive and negative A, 

p{2N, -A) = Up*{2N, A)U (5) 

p(2N + 1, -A) = T, v Up*{2N + 1, A)UH y , (6) 
where T, y — (cr 27 ) 18 ™, and U = [ ® cr^, and conjugate is taken in the basis where a z is diagonal. In view of the 

n odd 

above symmetries, below we restrict to the case of positive A > 0. 
We shall consider two particular limits: 

A) weak coupling, r<l, strong driving (large gradient) k = 1). 

B) weak gradient k<1, strong coupling <C 1 

Both limits show qualitatively similar singular behaviour as k — >■ or T — > 0, at the isotropic point A = 1 as 
demonstrated below. 
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III. CASE OF STRONG DRIVING k = 1 AND WEAK COUPLING Y < 1 



To investigate the system in the case of large boundary gradients n = 1 and weak coupling T <C 1, we solve the 
linear system of equations which determine the full steady state, analytically, making use of the global symmetry (U) 
[22j which decreases the number of unknown variables by roughly a factor of 2. For N = 2,3,4 we have indeed 9, 
31 and 135 real unknowns, respectively, instead of 2 2N — 1 real unknowns for the full Hilbert space. For N = 2 the 
current j(T, A) is readily obtained as 



i(r,A) 



8r 4 (3 + er 2 + 2A 4 ) 



24r 6 + 3(A 2 - l) 2 + 2r 4 (7A 4 + 30) + T 2 (2A 4 + 13A 2 + 30) ' 
from which we see that at the points A = ±1 and in the small coupling limit, it scales quadratically with T, 

m 



lim lim 



r^o A-Ti T 2 

Note, however, that for A ^ ±1 the current scales as T 4 , i.e 

m 



N=2 



lim 



N=2 



r^o r 4 

For N = 3 the current j(T, A) at A = 1, still scales as T 2 , 

m 



8(2A 2 + 3) 
3(A 2 - l) 2 



lim lim 

r-ioA-ii T 2 



N=3 



16 



(7) 



(8) 



(9) 



(10) 



but a different scaling behavior, with respect to N — 2 case, is obtained for other A points, where the j is found to 
scale also quadratically, 



lim 



r-K) r 2 

Taking the limit A — > 1 in the above, we obtain 



JV=3 



6A (9A 2 + 59) 
9A 4 + 155A 2 + 416' 



lim lim 

A->ir->o T 2 



N=3 



102 
145' 



different from (|10l) . Consequently, the limits A — > 1 and r — > do not commute. In more details, we have 



lim 
A->i r 2 



N=3 



16(3 + r 2 ) 

27 + 52r 2 + 12r 4 ' 



(11) 



(12) 



(13) 



ir=o 



from which Eq ([T0|) is straightforwardly obtained. Alternatively, one can see the presence of a singularity at all even 
orders of expansion of j(T): starting from the order 4 of the expansion, all even terms (d m jo(T, A)/dT r ' 
poles at A = ±1 of the type 1/(A 2 - l) m " 2 , e.g. 

d m Jo(r,A) _ = (-1) 

dT rn lr-0 ( A2 _ x y 



m+1 



r^Q(A) 



contain 



(14) 



with Q(A) a rational function which is regular at points A = ±1. 
For N = 4 one can see that the current at A = ±1 scales with T as: 



lim 



N=4 



8(27 + 46r 2 + ior 4 ) 

3(27 + 135r 2 + 144r 4 + 28r 6 ) 



a->±i r 2 

from which it follows that in the small coupling limit the current reduces to 

i(r) 8 



lim lim 

r-ioA^±i r 2 



(15) 



(16) 



N=4 
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At all other points A ^ ±1, however, the current scales as for the N = 2 case, e.g 



lim ™ 



r-s-o r 4 



N=A 



Although analytical expressions of the current are very difficult to obtain for larger values of N (e.g. A > 4), we 
can infer from numerical calculations that similar behavior exist also in these more complicated cases. Indeed, from 
the exact expressions reported above and from direct numerical calculations it is possible to extrapolate the values of 
the current for arbitrary N at the points A ± 1 in the small coupling limit as: 



lim lim 



r-»oA->±i r 2 



= §(#-!), (18) 

N y 



hinting at T = becoming a singular point in the thermodynamic limit. One can readily see that our conjecture (jT5J) 
coincides with the exact values reported in Eqs. flSJ), (fTU|) . (|TB)) for cases N = 2, 3, 4, respectively. Moreover, numerical 
results provides us a high confidence about the validity of Eq. (TT8"]) for arbitrary N. More precisely, we find from direct 
numerical solutions of the LME that the peaks of the current for N = 7,8, 9 at T = 1CT 5 are: 5.3333,6.2222,7.111, 
respectively, in perfect agreement with the prediction of Eq . (fT5|) . Also note from Fig. [T]that the numerical values 
of the current peaks depicted in panels (c),(d) for cases TV = 5,6 (3.5556 and 4.4444, respectively) are in excellent 
agreement with Eq. (TTB")) . Moreover, using the Matrix Product Ansatz [l5j], we were able to check the conjecture 
(p~8|) analytically up to sizes N = 100 (details to be published elsewhere). From the physical point of view, Eq. (fT8|) 
is quite interesting because it implies that in the weak coupling limit an addition of an extra spin to a finite chain 
contributes to the total current by a "quantum" of 8r 2 /9. Since the spin current is bounded, the Eq. tfTB")) hints at a 
singularity in T = point in the thermodynamic limit N — > oo. Further details and an analytical proof of (|18p will 
be presented elsewhere. 

One naturally expects that the noticed singular behaviour of the driven XX Z model at the isotropic point is not 
restricted to the magnetization current but can be seen for other observables as well. Here we show that the energy 
current is also anomalous at the isotropic point. Our purpose here, in addition, is to demonstrate, that the global 
symmetry @ is not crucial for singularity, because the energy current only flows if the symmetry ((3]) is broken. 
Indeed, the energy current operator J E (A), defined through the local Hamiltonian term h n , n +i as dh n , n +i/dt = 
J,f_ 1 (A)-J r f(A), 

Jn ( A ) = + A(j„_i j „CT^ +1 + a z n _^ n+1 ) 

where j„ jTO = 2(er^er^ — er^er^), changes sign under the action of (jl]), (RU^ ot T, x ) J E [RU^ ot T, x ) — —J® and therefore in 
the stationary state it is strictly zero ( Jjf ) = J E = (as a local integral of motion, the energy current does not depend 
on the choice of the bond n) . Lifting up the symmetric choice for the couplings to the left and right reservoirs in ([1} , 
i.e., considering the Lindblad equation (TTJ) with Y l ^ r# we break the symmetry ^ and the energy current becomes 
permissible, see also [2(|. The current of energy for small couplings r^r^ <C 1, with a fixed ratio Tl/Tr = const 
shows, again, an approach to a singularity in the isotropic point, but with a different scaling, J E (A = I) = O^Tl), 
and qualitatively different behaviour, see FigjSJ Indeed, the anomaly of energy current at the isotropic point has a 
shape of a twin peak, which, as Tl decreases, becomes more and more narrow, see Fig 2(c) and for — > develops 
a twin peak singularity. For N — 4 we see a twin peak singularity at A = 1 and a single peak singularity at A = 0.5, 
illustrated on Fig(^b). Unlike the spin current which is the odd (even) function of A for odd (even) number of sites, 
the energy current is always an odd function of A. These parity features are direct consequences of the mappings ([5]) 
and ([5]), see also Further details will be discussed elsewhere. Here we just mention another observable which has 
a singularity at the isotropic point, also in the symmetric case = Tr = V: it is the actual value of magnetizations 
at the boundaries: (erf ), (&%■}, (a\),(a v N ), which, as T — > 0, remain finite at A = I, and attain different values at 
the point Al ± 0, data not shown. Interestingly, this singularity cancels in the difference Ax = (erf — afj) and 
Ay = (a\ — u v N ): the actual magnetization differences Ax, Ay are regular at A = I. 

At this point one might suspect, that the singularities in the currents of spin and energy are just artefacts (or 
direct consequences) of the singular behaviour of the boundary gradients. To eliminate the influence of the irregular 
behaviour of the boundary gradients, in the next section we consider another limit, where we control not only the 
boundary gradients Ace, Ay, but also the individual magnetizations at the boundaries (af),{o~f), (af},(afj),(aff), 
{o~ z N ), which become exactly equal to 0, — k, 0, k, 0, 0, respectively. 
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FIG. 1: Typical dependence of renormalized current j/T 2 versus anisotropy for small couplings Y and odd system size (Panels 
(a),(c)) and for even system size (Panels (b),(d)). Dashed, thin, and thick lines correspond to T = 10 _1 , 10 -2 , 10 -3 respectively. 
The system size is N = 3, 4, 5, 6 for Panels (a),(b),(c) and (d), respectively. The curves are obtained by a numerical solution of 
the LME. 



IV. CASE OF WEAK BOUNDARY GRADIENTS AND STRONG COUPLING 

In this section we investigate the weak gradient and strong coupling case (k <C <C 1) by means of the 

perturbative approach for the LME in powers of developed in [19j . Note that while in the standard derivation of 
the Lindblad Master e quat ion one usually assumes the coupling to the reservoirs to be weakQ , 0] , in the ancilla Master 
equation construction |21| this restriction is lifted and the effective coupling can become arbitrarily strong. We also 
remark that the uniqueness of the nonequilibrium steady state for any coupling T is guaranteed by the com plet eness 
of the algebra, generated by the set of operators {H, V m , W m , V*, W^} under multiplication and addition [26|, and 
is verified straightforwardly along the lines [27j. We search for a stationary solution of the Lindblad equation in the 
form of a perturbative expansion in T -1 , 



p(A,r) = ^(-) Pk (A), (19) 



where Po(A) = limr_>.oo p{&, T) satisfies Clr[po(A)] = 0. Here and below we denote by £/,/{[.] the sum of Lindblad 
boundary actions = + £r[-]- This enforces a factorized form 



p (A, K ) = Pl ® r + A/ o( A - K ) ® PR, (2°) 




where and pa are one-site density matrices given by (O,© which satisfy Cl[pl] = 0, Cb{pr] = and Mq(A, k) is 
a matrix to be determined sclf-consistently later. Below we shall drop A- and iV-dependcncc in p k and Mk for brevity 
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FIG. 2: Dependence of renormalized energy current j/Tl versus anisotropy for small couplings r_L,r_R and system size N — 3 
(Panel (a)) and for JV = 4 (Panel (b)). Dashed, thin, and thick lines correspond to Tl — 0.1,0.05,0.02 respectively, while 
Tr = Tl/2. The curves are obtained by numerical solution of the LME. The peaks for JV = 4 are at the points A = ±0.5, ±1. 
Panel (c). A closeup view of the Panel(a), demonstrating double peak structure close to A = 1. 



of notations. Substituting (fl9|) into (fTJ), and comparing the orders of T fc , we obtain recurrence relations 



i[H, Pk] = ^LRpk+i k = 0, 1, 2... 



(21) 



A formal solution of the above is Pk+i = ^^\t{Qk+i) where Qk+i = —i[H,pk]- Note however that the operator 
Clr has a nonempty kernel subspace, and is not invertible on the elements from it. The kernel subspace ker(£ifl) 
consists of all matrices of type pl® A® pn where A is an arbitrary 2 N ~ 2 x 2 N ~ 2 matrix. Therefore Pk+i exists only 
if [H, pk] H ker(£i^) = 0, which in our case reduces to a requirement of a null partial trace, see [l9l |. 

Tr hN {[H,p k ]) = 0, k = 0,1, 2.... (22) 

which we call secular conditions. Finally, Pk+i is defined up to an arbitrary element from ker(££#), so we have 



Pk+i = 2£ L 1 R (i[H, p k }) +p L ® M k+ i ®p R , k = 0, 1, 2... 



(23) 



Eqs. (|20p . (|2"2")l and (f2"3")l define a perturbation theory for the Lindblad equation (TTJ) for strong couplings. At each 
order of the perturbation theory the secular conditions (|22|) must be satisfied, which impose constraints on Mk- Our 
aim here is to determine the exact NESS in the limit r — >• oo for which considering the two first orders k — 0, 1 has 
proved to be enough. The case N = 2 is trivial since Pk = PL® PR- For N = 3, the most general form of the matrices 
Mo(A), Mi(A), by virtue of the symmetry is Mo = bo(cr x — <r y ), Mi = bi(a x — a y ) where b$,bi are unknown 
constants. The secular conditions (|22l) for fc = do not give any constraints on bo, while those for fc = 1 give two 
nontrivial relations b\ = and — 4Ak 2 + (3 + 4A 2 — k 2 )6q = from which both Mq and Mi are determined. For 
= 4 each of the matrices Mo (A), Mi (A) compatible with the symmetry Q, contains 9 unknowns, all fixed by the 
secular conditions (|2"2")l for k = 0, 1, and so on. For JV = 3 wc obtain 



Mo (A) 1^ 



4Ak 2 
3 + 4A 2 - k 2 



((T X -<7 y ). 



(24) 



7 



The respectively stationary current ja(n) = limr-^oo j(^,r) 



Jo(k)U= 3 



8Ak 2 



3 + 4A 2 - k 2 



(25) 



does not have any non-analyticity as k — > 0, unlike in the case of weak coupling. However, for N > 3 the singularity 
sets in. For N = 4 the exact expression for the current is 



io(«) 



7V=4 



8A 2 k 4 (16A 4 + 2A 2 (5k 4 - 11k 2 + 24) - 8k 6 + 25k 4 - 36k 2 + 27) 

D(A,k) 



(26) 



where 



D(A, k) = 16A 8 (k 2 + 3) + 4A 6 (k 6 + k 4 + 48k 2 - 6) + A 4 (-6k 8 + 21k 6 - 
133k 4 + 339k 2 - 69) A 2 (9k 8 - 25k 6 - 7k 4 + 117k 2 + 18) 
-2k 8 - k 6 + 3k 4 - 27k 2 + 27, 

and it does contain a singularity at A = 1. Indeed, for N = 4 we observe a non-commutativity of the limits A — > 1 
and k — > 0, signalizing the presence of a singularity, 



lim lim JoiK) 







(27) 



JV=4 



lim lim 

t->0A4l K 
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(28) 



easily verifiable from (|26[) . Indeed, taking the lim K _j.o first, we expand j(k) in orders of k and find the first nonzero 
contribution at the fourth order, 



lim M 



*0 



8A 2 (4A 2 + 9) 

(29) 



N=4 



3(A 2 -l) z (4A 2 + 3) 



which is singular at A = 1. As a consequence, the current at the point A = 1 (and, in virtue of the symmetry ([5]). 
also A = —1) has a different scaling (jo(K, A = ±1) ~ k 2 ) then all other points A ^ ±1, where the current scales as 
k 4 , see Fig|3l Note that the singularity type for weak driving k is exactly the same as the one described in Sec lIIII for 
even number of sites and weak coupling T. 

For odd-sized system N = 5 the exact expression for the current is very complicated and some limiting cases are 
reported in the[AJ Analyzing the analytic expression in the limit k < 1 we find the scaling of the type jpp =5 (n) = 
q:(A)k 2 where the prefactor is a(A) is singular at A = 1, namely a(l) — > 8/7 and a(l ± 0) — > 64/181 as k — > 0, see 
also FigJ3J This type of singularity is exactly the same as the one seen in Sec lIIII for vanishing T and odd system sizes. 

As the system size increases, the order of polynomials grow with system size and becomes complicated. However, 
at the qualitative level we see the same behaviour as discussed above for N = 4 ( N = 5 ) for system of even (odd) 
sizes. 

The qualitative similarity in the magnetization current behaviour in the cases (i) of weak coupling and large driving, 
and (ii) weak driving and large coupling is not a trivial one. In the case (ii), the amplitude of the effective x- and 
y- boundary gradients is exactly equal to k, independently on the anisotropy A, so in the limit k <C 1 it becomes 
infinitcsimally small. Also the individual boundary magnetizations (a") ,(a^f) ,a = x,y,z are bounded by k. 

In the case of weak coupling r <SC 1, the amplitude of the boundary gradient at the isotropic point also vanishes in the 
limit r — > 0, but the individual boundary magnetizations (erf ),((?%) for a = x,y remain finite and are discontinuous 
at A = 1. 

The presence of the singularity and non-commutativity of the limits limA->i and limr-!.o(or lim K _s.o) must be 
related to the existence of an additional symmetry which the XX Z model acquires at the isotropic point A = 1. This 
symmetry however is not the full rotational symmetry of a unitary quantum XXX Hamiltonian, since the dissipative 
Lindblad terms in LME do not have the full rotational symmetry On the other hand, qualitative differences between 
singularities for odd and even system sizes TV, seen for spin current, energy current, and other observables, is a 
consequence of intrinsic properties of the XY-twisted model, which is manifested by ([SJ) , (JSJ) - 
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FIG. 3: Typical dependence of renormalized current j'/k 2 versus anisotropy for small driving k and even system size (Panels 
(a)) and for odd system size, Panel (b). Dashed, thin, and thick lines correspond to k = 0.2, 0.1, 10~ 2 respectively. The system 
size N = 4 for Panels (a) and N — 5 for Panels (b). The curves are given by exact analytic expressions (|26p for N = 4 and 
analogous analytic expression (not reported) for N = 5. 



V. CONCLUSIONS 



We have investigated an XX Z spin chain coupled at the ends to a dissipative boundary reservoirs, which impose 
a twisting angle between the first and the last spin in the XY plane, described by a Lindblad Master equation. We 
pointed out the non-analytic character of the non-equilibrium steady state in the two limits: of vanishing coupling 
r and of vanishing driving k. Unlike in the approaches proposed before, where a small but fixed boundary driving 
(of the order of 10 -2 ) is being typically used, here we access arbitrary small coupling or driving analytically. Our 
approach allows to establish a presence of a singularity in the NESS for T — > and n — > 0, which is expected to be 
present in the system for all system sizes. The singularity is evidenced on an example of several observables: the 
magnetization current, the energy current, and the boundary magnetizations. For magnetization current, the analytic 
treatment is presented. 

The character of the singularity qualitatively depends on a parity of system size N. For odd N we find that the 
spin current scales for small T or small k as j — /(A)T 2 , j = g(A)n 2 , where for vanishing T or k the functions / and 
g have different finite values at A = 1, and A = 1 ± 0. For even sizes N, the spin current still scales quadratically 
with r or k at the point A = 1, while it scales as j = fi(A)T i ,j = <7i(A)k 4 at all other points. The energy current 
J E scales in the isotropic point linearly with couplings r^r^ (which must be different, Tl ^ ^r), and develops a 
twin peak singularity. Here we are not already in a position to discuss the conductivity of spin or energy due to the 
small system sizes; however the DMRG studies [HI and recently proposed exact approaches might give access to large 
system sizes and even to the thermodynamic limit. 
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through a FIRB initiative. M.S. acknowledges support from the Ministcro dell' Istruzione, dell' Universita c dclla 
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Appendix A: Current for N = 5 (weak driving limit) 

For N = 5 and in the limit r — > 00 the exact expression for the current is given by a complicated expression. We 
report different limits here below. Taking the limit lim K _;.o first, we find 

ic^A) _ 16A(A 2 + 3) 
k-To K 2 28A 4 + 81A 2 + 72' 1 1 

and consequently limA->i lim K _yo J "^ A ^ = On the other hand, taking the limit limA->i first, we find 

limA-»i J "^ A ^ = — 77, where 
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F = 8(128k 22 - 2128k 20 + 57 92k 18 + 44096k 16 - 245827k 14 + 182930k 12 + 
980805k 10 - 1523750k 8 - 1515697k 6 + 3477286k 4 - 1171473k 2 + 2235870), 



G = 576k 22 + 856k 20 - 34736k 



lcS 



102787k 16 + 1669223k 14 - 4033065k 12 - 



2004745k 10 + 13257113k 8 + 3528813k 6 - 16075891k 4 - 7450779k 2 - 15651090 




liniA-s-i J °^ A ' 1 = |. Another way to see the discontinuity is to estimate the derivatives 
all even orders of which, starting from m = 4, have singularities at A = 1 as 




For odd m = 2n + 1, we find (<9 2n+1 j (K, A) /Ok 



= 0. 
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